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Abstract. We study the spherical collapse model in dark energy cosmologies, in which dark energy is modelled as a mini- 
mally coupled scalar field. We first follow the standard assumption that dark energy does not cluster on the scales of interest. 
Investigating four different popular potentials in detail, we show that the predictions of the spherical collapse model depend 
on the potential used. We also investigate the dependence on the initial conditions. Secondly, we investigate in how far pertur- 
bations in the quintessence field affect the predictions of the spherical collapse model. In doing so, we assume that the field 
collapses along with the dark matter. Although the field is still subdominant at the time of virialisation, the predictions are 
different from the case of a homogeneous dark energy component. This will in particular be true if the field is non-minimally 
coupled. We conclude that a better understanding of the evolution of dark energy in the highly non-linear regime is needed in 
order to make predictions using the spherical collapse model in models with dark energy. 
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Q ' General Relativity, needs to have negative pressure. 



1. Introduction 

The most surprising result of observational cosmology is the discovery that high redshift supernovae are less bright than expected 
[Perlmutter et al (1999)| [Tonry et al. (2003), The interpretation of this finding is that the expansion of the universe accelerates, 
instead of slowing down. The matter responsible for this accelerated expansion is dubbed "dark energy" and, according to 
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One important goal of cosmology is to obtain information about the nature of dark energy (see e.g. |Ratra & Peebles (2002)] 
for a review). The physical origin of dark energy is unknown and its likely that new physics beyond the standard model of particle 
physics will be needed to explain its properties. Scalar fields are plausible candidates for dark energy (see e.g. |Wetterich (1988a)| 
IWetterich (1988b)||R'atra & Peebles (1988)^.Wetterich (1995)^ (these scalar fields are called quintessence fields). A distinctive 
feature of these models is that the properties of dark energy, such as its equation of state, generally vary during the cosmic 
history. 

Cosmological observations probe the properties of dark energy. The anisotropy spectrum of the Cosmic Microwave 
Background (CMB), the shape of matter power spectrum and the distance-redshift relation are possible sources of information 
and will help us to distinguish between the different models of dark energy. 

The behaviour of small perturbations in a scalar field and its effect on CMB anisotropics and structure for- 
mation has been investigated by a number of authors. However, the behaviour of quintessence during the grav- 
itational collapse into the highly non-linear regime is not well understood and currently under investigation 
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Usually it is assumed that there are no density fluctuations in the quintessence field on cluster scales and below. The reason 
for this assumption is that, according to linear perturbation theory, the mass of the field is very small (the associated wavelength 
of the particle is of the order of the Hubble radius) and, hence, it does not feel overdensities of the size of tenth of Mpc or smaller 
[Wang & Steinhardt (1998)^ Let us assume for instance the particle physics candidate for dark energy, a scalar field <p. In the linear 
regime of the small cosmological perturbations, it can be shown that, during the matter dominated epoch, perturbations in the 
scalar field are described by [Hwang & Noh (200 1)| 

6if> + 3H64> + (k^/a^ + V")6(f> = 4>Srdm (1) 
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where the metric was perturbed about a flat Friedmann-Robertson- Walker (FRW) metric. From equation Q the effective Jeans 
length for linear perturbation of the scalar field is roughly given by 

Aj ^ In I y/V^^ (2) 

which turns out to be of the order of the horizon size |Ma et al. (1999)| Thus, on scales much smaller than Aj the fluctuations 
on the field are unimportant. This results lead people to assume that dark energy is an homogeneous field throughout the uni- 
verse, which is not affected by structure formation on small scales. However, the source term 4>Scdm implies that dark energy 
is not smoothly spread [Caldwell et al. (1998)| And in spite of the perturbations in dark energy are present and grow only on 
scales about the horizon size and larger, the effects on the evolution of matter overdensities are indeed important and significant 
Ferreira & Joyce (1997) For instance, the spatially inhomogeneous distinctive property of quintessence, with respect to the cos- 
mological constant A, has an important effect on the large-angular-scale of the CMBR [Steinhardt (2003)| To study the effect 
of dark energy on small scales, N-body simulations were used and the assumptions made that quintessence perturbations are 
negligible. (In Maccio et al. (2003b), a possible coupling of dark energy and dark matter was explicitly taken into account for the 
perturbations.) 

The assumption of neglecting the effects of matter perturbations on the evolution of dark energy (and its backreaction) at 
small scales is indeed a good approximation when perturbations in the metric are very small. Notice however, if the field couples 
explicitly to matter, such as in the coupled quintessence scenario (see e.g. |Wetterich (1995)| and |Amendola (2000)t , the conclusion 
that the Jeans length is of order the horizon size should no longer hold. 

Nevertheless, even if the field is not coupled, one should be careful when extrapolating the small-scale linear-regime results 
to the highly non-linear regime. Then, locally the flat FRW metric is not a good approximation anymore to describe the geometry 
of overdense regions. It is natural to think that once a dark matter overdensity decouples from the background expansion and 
collapses, the field inside the cluster feels the gravitational potential inside the overdensity and its evolution will be different from 
the background evolution. The backreaction effects in the highly non-linear regime could influence the evolution of perturbations 
in dark energy considerably, which in turn influence the evolution of the matter perturbation. In |Bean & Magueijo (2002)| it was 
suggested that the quintessence field could have an important impact in the highly non-linear regime. Due to these considerations, 
one may even ask if the quintessence field can be important even on galactic scales? It was speculated by [Wet terich (20 1)| 
IWetterich (2002)1 and by |Arbey et al (200 1)| that highly non-linear perturbations might indeed be important even on galactic 
scales. It was found that, at least in principle, the quintessence (or a scalar field) field could be responsible for the observed flat 
rotation curves in galaxies. In Padmanabhan & Cho udhury (2002)] [Padmanabhan (2002)| [Bagla et al (2003)| and |Causse (2003)| 
more exotic models, based on tachyon fields, have been discussed and it was argued that the equation of state is scale-dependent. 

If it turns out that backreaction effects of metric and density perturbations in dark matter could influence perturbations of 
quintessence on small scales, this could significantly change our understanding of structure formation on galactic and cluster 
scales in models with quintessence. 

One popular model to study the non-linear growth in cold dark matter is the spherical collapse model (see. e.g. 
[Padmanabhan (1995 )| or [ Peacock (1999)) . The model was used first in the standard cold dark matter scenario, but later also 
in the cold dark matter model with cosmological constant (ACDM) by |Lahav et al. (1991)| Recently, the model was also used 
in quintessential scenarios and dark matter under the assumption that dark energy does not cluster on scales much smaller than 
the horizon (see e.g. [Wang & Steinhardt (1998)). It was then subsequently use d in order to make predictions for cosmological 
observations (see e.g' |Weller et al. (200 1)| [Weinberg & Kamionkowski (2003)|[Battye & Weller (2003)) . 

The aim of this paper is two-fold. Firstly, we model dark energy as a non-minimally coupled scalar field and consider four 
different potentials in order to investigate how the predictions of the spherical collapse model depend on the potential. In addition, 
the effect of the initial conditions are investigated as well. Our second aim is to get a feel about the effects of perturbations in the 
quintessence field on the predictions of the spherical collapse model. In doing so, we investigate two extreme cases: in the first 
case we follow the literature and assume that there are no fluctuations in the quintessence field on the scales of interest. We then 
assume the other extreme case and assume that the field inside the overdensity collapses together with the dark matter' . As we 
will see, the predictions of the spherical collapse model will be quite different in this case, although the field is only slightly non- 
linear (with density contrasts of order one) at the time of virialisation, whereas dark matter is in the highly non-linear regime. Our 
results imply that the predictions of the spherical collapse model depend on the assumptions made for the clustering properties of 
the quintessence in that model. Thus, a better understanding of the behaviour of the scalar field in the highly non-linear regime is 
needed, if the spherical collapse model is used for predictions involving large scale structures, such as cluster abundances, weak 
and strong lensing, etc. 

The paper is organised as follows: In Section 2 we describe briefly the spherical collapse model and write down the equations 
used. We also give an overview of the potentials used. In Section 3 we describe our assumptions and give the results of our 
numerical calculations. A discussion of the results and our conclusions can be found in Section 4. 



' This approach was first used, in the context of varying-fine structure theories [Vlota & Barrow (2003)||Mota & Barrow (2004)| 
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2. The spherical collapse model 

We consider a flat, homogeneous and isotropic background universe with scale factor a(t). Since we are interested in the matter 
dominated epoch, when structure formation starts, we assume that the universe is filled with cold dark matter of density p„i oc a^^ 
and a dark energy fluid, with energy-density p^. The equations that describe our background universe are (we set h - c = 1 
throughout the paper): 

3//2 = SnG{p,n+p^) (3) 
= -3H(l + w^)p^ (4) 



where H = a/a is the Hubble rate. When - -1 then dark energy is the vacuum energy density A. If dark energy is a scalar 

+ y(0) and = \cf 



field <p (quintessence), p^ = -i- V((p) and = - V((f>), where V((f>) is the scalar field potential. In this case, it is useful to 



re-write equation Q as 

ip + 3H^ + y' = (5) 
where the prime represents a derivative with respect to (p. In this paper we consider four examples of quintessential potentials: 

- The double exponential potential fBarreiro et al. (2000)| 

y(0) = M (exp(;80) + exp(r0)) . (6) 

- The exponential potential with inverse power |Steinhardt et al. (1999)| 

y(0) = M(exp(r/0)-l) (7) 

- The Albrecht-Skordis model |Aibrecht & Skordis (2000)| 

V(cf,) = M (a + (0 - B)2) exp(-r0) (8) 

- The supergravity-motivated potential ^Brax & Martin (1999)| 

y(0) = Mexp(02)/0>' (9) 

We choose the parameters in the potentials and the initial conditions in the background such that today we obtain Q.„, - 0.3, 
= 0.7, Hq = 100/ikm/(Mpc- sec) with h = 0.7 and -1 < < -0.8. 

In order to study the non-linear evolution of the gravitational collapse we make use of the spherical collapse (or infall) model. 
Here, one considers a spherical overdense region of radius R and models the interior spacetime as a FRW universe. This approach 
is equivalent to study the effect of perturbations to the Friedmann metric by considering spherically symmetric regions of different 
spatial curvature in accord with Birkhoff 's theorem. Clearly, this model ignores any anisotropic effects of gravitational instability 
or collapse. However, the model is quite useful to learn about non-linear gravitational collapse and it is widely used as a starting 
point for (semi-) analytical models of large-scale structure formation. 

Consider a spherical perturbation, in the cold dark matter fluid, with (spatially) constant internal density which, at an initial 
time, has an overdense amplitude 5, > and «: 1. The cold dark matter density inside the cluster is then initially pcdm = 
Pm(l + S(Ri, ti)), where Rj is the initial radius of the overdensity. The dark energy fluid density inside the overdensity, p^ , will 
initially be considered the same as the background one^. At early times the sphere expands along with the background, but with 
a slightly different expansion rate. For a sufficiently large 5,, gravity prevents the sphere from expanding for ever: at one point, 
the overdensity will stop expanding and start to collapse. Three characteristic phases during the evolution can then be identified: 

- Turnaround: the sphere breaks away from the general expansion and reaches a maximum radius. 

- Collapse: if only gravity is significant, the sphere will then collapse towards a central singularity where the densities of the 
matter fields would formally go to infinity. In practise, pressure and dissipative physics intervene well before this singularity 
is reached and convert the kinetic energy of collapse into random motions. 

- Virialisation: dynamical equilibrium is reached and the system becomes stationary: the radius of the system the energy of 
the different components are constant. 

In the spherical collapse model, due to its symmetry, the only independent coordinates are the radius of the overdensity and 
time. Also, as is standard practise when using this model, we consider that there are no shell-crossing which implies mass (of cold 
dark matter) conservation inside the overdensity and independence of the radius coordinate jPadmanabhan (1995)| The equations 
can then be written ignoring the spatial dependence of the fields (but still including an equation for the evolution of the radius). 
The evolution of a spherical overdense patch of scale radius ^(f) is given by the Raychaudhuri equation: 

3R = -AnGR (pcdm + P0,(1 + 3w^^) (10) 
^ This assumption will be made throughout the paper. 
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Note that it would be a wrong to use the Friedmann equation for a closed universe with a constant curvature k, since the former 
can vary in time [Weinberg & Kamionkowski (2003), Wang & Steinhardt (1998). In the cluster, the evolution of p^^ and pcdm is 
given by 

= -3^(i + w^jP0, + r (11) 

and Pcdm °^ due to mass conservation. The quantity F describes the energy loss of dark energy inside the dark matter 
overdensity, as this component does not necessary follow the collapse of dark matter and energy can formally flow out of the 
system. As such, F encodes in how far backreaction effects from the dark matter non-linearities act on the dark energy component. 
Formally, F can also describe the coupling between matter and dark energy, in which case we would have to add a term containing 
F in the energy conservation equation for dark matter In this case it would be no longer true that pcdm °^ R^^ ■ 
Once again, in the case of a scalar field, equation (II U can be written as 

^, + 3^(j)c + Vc'{<pc)^l; (12) 
K <p 

where (pc is the field inside the overdensity in order to distinguish it from the background value and Vc = V(0e) is its potential. 

We will evolve the spherical overdensity from high redshift until its virialisation occurs. According to the virial theorem, 
equilibrium will be reached when T - jRj^', T - jMv^.^ is the total kinetic energy at virialisation and v^^.^ is the mean-square 
velocity of the components of the cluster, and U is the average total potential energy in the sphere. It is useful to write the 
condition for virialisation to occur in terms of the potential energies associated the different components of the overdensity. The 
potential energy for a given component 'x' can be calculated from its general form in a spherical region |Landau & Lifshitz (1975)| 

t/, = 27T p,o,^,r^dr, (13) 



(r) = -27rGp,(l + 3w,) - y j , (14) 

where p^, is the total energy density inside the sphere, O ^ is the gravitational potential due to component x with energy density 
p t and equation of state Wx- In our calculation, we have taken also the pressure contribution from dark energy into account, since 
the Poisson-equation reads (including the correction from General Relativity) 

AO = 47rG (p + 3f ) . (15) 

Using the virial theorem and energy conservation at the turnaround and cluster virialisation times, we obtain an equilibrium 
condition only in terms of the potential energies: 

1 d 

■fo^'^Uc + f/*) + Uc + t/* lz„ = f/c + (16) 

where Uc is the potential energy for Cold Dark Matter, U^^. is the potential energy for dark energy, z,, is the redshift of virialisation 
and Zta is the redshift at the turnaround of the over-density at its maximum radius, when R - R,„ax and = 0. In the case 
where = -1 the expressions reduce to the usual virialisation condition for ACDM models. One should point out here the 
inconsistency widely used in the literature, when one makes use of equation M6\ together with the assumption that energy is not 
conserved inside the overdensity (dark energy is assumed to be homogeneous, so it has to flow out from the overdensity). In the 
case where F = the whole process of virialisation is self-consistent. 

In addition, we also discuss the linear growth factor D{z) - {Sdz) / a{z)) / (6^0) / a{Q)) for the four different potentials. The 
linear density contrast 6c obeys the equation 

6c + 2H6c - 4-nGp„,6c =0. (17) 

As we will see, the growth factor D(z) depends on the potential used. 

The behaviour of dark energy during the evolution of a cluster can now be obtained by numerically evolving the background 
eqs.Q-llSli and the cluster eas.( ll0t -( ll2t until the virialisation condition il6\ holds. Additionally to the four potentials mentioned 
above, we also have considered three different types of large scale structure models: the Standard Cold Dark Matter model 
SCDM, the ACDM model, and a model with dark energy with constant equation of state - -0.8 (QCDM). The SCDM and 
the ACDM model was used to test our code and we will not give the results for these models, as they are well known. 

In next section we will make two assumptions about the function F. In one case we assume that the field is homogeneous all 
over the Universe, so that the field inside the overdensity evolves just like the background. In the second case we assume that the 
field collapses together with the dark matter. We would like to point out that our assumptions are different from the ones made 
in the first version of'Loka s & Hoffman (200 1}| where it was assumed that dark energy is homogeneous but a Friedmann-like 
equation was used to study the evolution of Rit). 



David F. Mota and Carsten van de Bruck: On the spherical collapse model in dark energy cosmologies 



5 



3. Non-linear collapse with dark energy 

We now study numerically the equations given in the last section. For this, we must make assumptions about the function F, 
which encodes the backreaction of dark matter perturbations on dark energy. We will make two assumptions for F: in the first 
case, the quintessence field is assumed to be smooth throughout space. This can be obtained by putting 

r = -3(---)^? (18) 

and 0c(fi) - (piti) ™d (pdti) - (piM)- Then we obtain 0c(O - <l>{t) for all times. This is the standard assumption made in the 
literature. As already mentioned in the introduction, the main reason for this assumption is that the effective Jeans length for 
linear perturbations turns out to be of the order of the horizon size. Thus, so the argument goes, on scales much smaller than /Ij 
the fluctuations are unimportant. 

The form of F chosen above is clearly not the full truth. At least at very late times during the collapse of the dark matter, 
especially when the density contrast in dark matter is very large (<5m » 1), the field should no longer feel the background metric, 
i.e. expand with the background, but decouple from it. In this regime, the evolution of dark energy could be different and influence 
the details of the collapse. The detail of this can only be obtained from a fully relativistic calculation, which is beyond the scope 
of this paper. However, to get a feel of what may happen we make now the other extreme assumption that the field follows the 
dark matter collapse from the very beginning. That is, we assume that F = 0. This is clearly not a realistic assumption either, but 
the results we obtain are surprising and interesting in itself. 



3.1. Results for the four quintessence potentials 

We begin our discussion with an analysis of the four different quintessential potentials mentioned in Section 2. The results of our 
numerical calculations are presented in Fig. [2 13 and |5 




Fig. 1. Quintessence model V = M(exp(/3(l>) + exp(-ytf>)) \Barreiro et al. (2000)\ Top left panel: Evolution of in the background (dashed 
line) and inside an overdensity (solid line) for T = Q, as a function o/log(l + z) (overdensity virialises at z = 0). Top right panel: Evolution 
P<l>/PiPc (solid line) and pmlPcdm (dashed line) in the case T = Q as a function o/log(l + z) (overdensity virialises at z = 0). Middle left 
panel: The ratio R,a/Rv in the case of an inhomogeneous (T = Oj (solid line) and homogeneous scalar field (dashed line). The ratio Pip/pmatta 
inside the overdensity (solid line j f F = 0) and in the background ( dashed line ). Bottom left: as a function of Zv, considering the effect of 
an inhomogeneous quintessence field (T = Q) (solid line). The results for a homogeneous quintessence field are shown as well (dashed line). 
Bottom right: Predictions for the linear growth factor for this potential (solid line). The dashed line represents the ACDM case. 
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Fig. 2. Quintessence model V = M(exp(y/c/>) - l) \Steinhardt et al. (1999)\ Top left panel: Evolution ofw^ in the background (dashed line) 
and inside the overdensity (solid line) for T = Q as a function o/log(l + z) (overdensity virialises at z = Oj. Top right panel: Evolution 
Pit>lp<Pc (^ol'd line) and p„/pcdin (dashed line) in the case T = as a function o/log(l + z) (overdensity virialises at z = 0). Middle left 
panel: The ratio R,a/Rv in the case of an inhomogeneous (T = Q) (solid line) and homogeneous scalar field (dashed line). The ratio Pift/pmnttet 
inside the overdensity ( solid line ) (T - 0) and in the background ( dashed line ). Bottom left: as a function of zv, considering the effect of 
an inhomogeneous quintessence field (T = 0) (solid line). The results for a homogeneous quintessence field are shown as well (dashed line). 
Bottom right: Predictions for the linear growth factor for this potential (solid line). The dashed line represents the ACDM case. 



In these plots we show for all models the behaviour of the equation of state inside the cluster and in the background (for the 
case r = 0), the density contrast in matter and field (again for F - 0), the ratio of the radius at turnaround Rta and radius at 
virialisation Ry for both the case F = and an homogeneous scalar field, the ratio of the energy densities in field and matter, the 
non-linear density contrast for both cases and the growth factor for the models compared to the ACDM model. 

As it is clear from the upper right plots, even in the F = case (ie. the quintessence field collapses with together with matter), 
the field will be subdominant in the virialised object. This can also be seen from the Figures on the right hand side in the middle 
panels in Figures m2]|3] and |4]: in all cases the ratio p^/pcdm is much smaller than one. Well before turnaround, i.e. while the 
density contrast in dark matter is smaller than one, our results agree with what one would obtain from linear perturbation theory. 
For F = 0, the perturbations in dark energy remain very small up to this point. Even after turnaround, our results imply that 
fluctuations in the dark energy component remain small, while dark matter enters the highly non-linear regime. 

Concerning the clustering properties of the quintessence field, it can significantly alter the predictions for the non-linear 
density contrast Ac = Pcdm,ciu,stei(fv)/Pcdm,backgi(fv) at virialisation: assuming that the field is inhomogeneous F = 0, for small 
virialisation redshifts the predictions can differ by a factor of four or more. In the same way the predictions for the ratio Rta/Ry 
strongly depends on how the scalar field behaves during the highly non-linear regime. The reason for this is, that the just 
before the overdensity formaUy becomes a singularity, the quintessence field becomes more and more important and changes the 
evolution of R(t) during the last stages of the collapse. We remind the reader that, in the Einstein-de Sitter model Ac ~ 147^ At 
high virialisation redshifts, all models predict Ac ~ 147, only at low virialisation redshifts, significant deviations can be expected 
due to the fact that the dark energy becomes more and more important. Note that in the Albrecht-Skordis model, the field behaves 
like a cosmological constant in the background (see FigO- Hence, one would expect that the differences to the ACDM model 
are small. This is indeed the case, because 0c ~ and, hence, F 0, as it can be seen from eq. ( I18> . Thus, in this model, the 
fluctuations in the quintessence field remain small and the field is almost homogeneous. 

^ If we would define A;.2 = Pcdm,ciuster('v)/Pcdm,backgr('c)> where ?c is the collapse time, its value in the Einstein-de Sitter model would be 
approximately 178. However, for computational convenience, we consider instead the quantity Ac at the time of virialisation. Ac and Ac2 differ 
only by a factor from each other. 
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Fig. 3. Quintessence model V = M(A + ((fi - B)^) exp(-'y^) [Albrecht &■ Skordis (2000)\ Top left panel: Evolution of in the background 
(dashed line) and inside the overdensity (solid line) for T = Q as a function o/log(l + z) (overdensity virialises at z = OJ. Top right panel: 
Evolution of p^lp^^. (solid line) and p„Jpcdm (dashed line) in the case T = Q as a function o/log(l + z) (overdensity virialises at z = OJ. 
Middle left panel: The ratio R,a/Rv in the case of an inhomogeneous (T = 0) (solid line) and homogeneous scalar field (dashed line). The ratio 
piplpmzua inside the overdensity (solid line) (T = Oj and in the background (dashed line). Bottom left: as a function of zv: considering the 
effect of an inhomogeneous quintessence field (T = 0) (solid line ). The results for a homogeneous quintessence field are shown as well ( dashed 
line). Bottom right: Predictions for the linear growth factor fi>r this potential (solid line). The dashed line represents the AC DM case. 



In Figure [S] we have plotted the case for a dark energy model with constant equation of state w = -0.8. As it 
can clearly be seen, the predictions for and Ry/Ria for the cases of a homogeneous and an inhomogeneous dark en- 
ergy component are very small in this case. This implies, that for models with constant equation of state and as long 
as the equation of state doesn't differ too much from w = -1, the fitting formulae presented in the literature are valid 
IWang & Steinhardt (1998)1 [Weinberg & Kamionkowski (2003)| 

The reason for the differences observed is that in the highly non-linear regime the field will play an important role in 
determine the collapse time if F = 0. Although at the time of virialisation the field is only subdominant, it will play an important 
role just before collapse. Thus, backreaction effects could significantly alter the predictions of the spherical collapse model. 

3.2. Dependence on the Initial Conditions 

We turn now our attention to the question, if the results of the spherical model depend on the initial conditions of the 
field (p (a discussion of initial conditions in the linear regime can be found in .Abramo & Finelli (2001). Doran et al (2003)| 
IBartolo et al (2003)t . We have run our code for different initial conditions and will discuss the case for an inhomogeneous 
(F - 0) and homogeneous dark energy component independently. The case for a homogeneous dark energy component is shown 
in Figure|6l In the left upper panel we show the evolution of the equation of state in the background for different initial values for 
at a redshift of lO"*. As it can be seen, the ratio p^/pcdm.cius is almost independent of the initial conditions chosen. The same is 
true for the predictions of Ac and Rt^/Ry. 

The case for F = is shown in Figure We find that if quintessence becomes inhomogeneous together with dark matter in 
the strong non-linear regime, the predictions of the spherical collapse model depend generically on the initial conditions of the 
field. 

Different initial conditions lead to different time evolution of the equation of state and the energy density of quintessence, 
both in the background and inside the overdensity (depending on the value of F). These two dark energy properties will affect the 
dynamics of the overdensity through equations ilOi and il4\ . hence affecting the predictions for the non-linear density contrast 
and the virialisation radius of the overdensity. 
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Fig. 4. Quintessence model V = M exp(0^ \Brax & Martin ( 1 999 j] Top left panel: Evolution ofw^ in the background (dashed line) and 
inside the overdensity (solid line) for T = Q as a function o/log(l + z) (overdensity virialises at z = 0). Top right panel: Evolution ofp^/p^^ 
(solid line) and p,„IPcdm (dashed line) in the case T = Q as a function o/log(l + z) (overdensity virialises at z = 0)- Middle left panel: The 
ratio R,a/Rv in the case of an inhomogeneous (T = 0) (solid line) and homogeneous scalar field (dashed line). The ratio p^lpm-ma inside 
the overdensity (solid line) (T = Oj and in the background (dashed line). Bottom left: as a function of zv: considering the effect of an 
inhomogeneous quintessence field (T = Oj (solid line). The results for a homogeneous quintessence field are shown as well (dashed line). 
Bottom right: Predictions for the linear growth factor for this potential (solid line). The dashed line represents the ACDM case. 



3.3. Comparison among the different models 

What are the differences in the predictions for the four quintessence models? In Figure|8]we plot the predictions for Ria/Rv and 
Ac for the different potentials with the different assumptions for F. As it can be seen from that Figure, the differences can be quite 
large if F = 0. If the field is homogeneous throughout space, the differences are small. Again, the reason is that the dynamics of 
the overdensity (see equations ilO\ and il4\ ) strongly depends on the dark energy properties (Potential, F, etc). 

In Figure|9]we plot the predictions of the linear growth factor {5c I a), normalised today, for the different potentials. As it can 
be seen, even in the linear regime there are differences between the four potentials. Such a behaviour was akeady noticed in 
IMaccio et al. (2003a)| 

3.4. Delayed collapse of dark energy 

The results presented imply that if dark energy collapses together with dark matter, it has an important impact on the predictions 
of the spherical collapse model. Clearly, our assumption on the parameter F so far is unrealistic. To gain more insight in how 
far our choices of F changes the predictions of the spherical collapse model, we make now the assumption that the field inside 
the cluster is the same as in the background until the turnaround, from which on it collapses together with the dark matter. 
We will call this situation as the delayed collapse of dark energy. In the light of results obtained in the linear theory (see e.g 
[Ferreira & Jo yce (1997)^ , this is a more realistic situation than studied up to now. The results for the precictions of A^ are shown 
inFig.Hni 

In Fig.[roiwe show the results for A^ for all four potentials. The solid lines are the cases for delayed collapse of dark energy, 
whereas the dashed lines correspond to the case of full collapse. As it can be clearly seen, the predictions for A^ change. Note, 
however, that the predictions are still far off from the the results obtained in the case of a homogeneous scalar field. Of course, 
the field is in the hnear regime, but it has still an important impact on the predictions of the spherical collapse model. 



David F. Mota and Carsten van de Bruck: On the spherical collapse model in dark energy cosmologies 



9 




log ( 1 + zj 




log ( 1 + z J 



E 



\ 

\ 

\ 

- \ 

\ 

\ 


^ ~ — — — ^ 




1 


1 


0.5 1 


1.5 


2 2.5 

log ( 1 + z ) 


3 


3.5 



Fig. 5. Model with constant equation of state w = —0.8. The upper panel shows the predictions for t^c for the inhomogeneous case (dashed 
line } and homogeneous case (solid line ). The middle panel shows the predictions for Ry /i?ta for the inhomogeneous case ( dashed line ) and 
homogeneous case (solid line). The lower panel shows the evolution of the density contrast in matter (dashed line) and dark energy (solid line) 
for the case of an inhomogeneous dark energy component (cluster virialises at z = 0). 
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Fig. 6. Dependence on the initial conditions for the double exponential potential \Barreiro et al. (2000)\ for the case of homogeneous dark 
energy. The solid and dashed lines show the results for different initial conditions. Top left panel: Evolution of in the background as a 
function o/log(l + z). Top right panel: The radio p^/pmatter inside the overdensity as a function of log(l + Zv). Bottom left: Evolution of os a 
function o/log(l + z,,)- Bottom right: Evolution ofR,„/Ry as a function o/log(l + Zv). Whereas the evolution of the equation of state is different, 
the predictions for I p^i^, and R^jR^ are only weakly dependent on the initial conditions. 
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Fig. 7. Dependence on the initial conditions for the double exponential potential \Barreiro et al. (2000)\ in the case for an inhomogeneous dark 
energy component (T = Q). The solid and dashed lines show the results for different initial conditions. Top left panel: Evolution ofw^ in the 
overdensity as a function o/log(l + z). Top right panel: The radio p^/pmaitei inside the overdensity as a function o/log(l + z,.). Bottom left: 
Evolution of Ac as a function o/log(l + Zv)- Bottom right: Evolution of R,a/Rv as a function o/log(l + Zv)- It can clearly be seen that the 
predictions in this case depend on the initial conditions. 





Fig. 8. The predictions for Rf^/R^ and for the different qiuntessential potentials. The left panels show the case for an inhomogeneous 
dark energy component (T - 0), whereas the right panel shows the results for the case of an homogeneous dark energy component. It can 
be seen that the predictions depend on the potential used, although the dependence is weaker in the case of an homogeneous dark energy 
component. Solid line \Barreiro et al. (2000}\ dashed-line \Steinhardt et al. ( 1999)\ dashed-dotted line [Albrecht &■ Skordis (2000)\ dotted-line 
\Brax &■ Martin (1999)\ 
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Fig. 9. The linear growth factor for the different quintessential potentials. It can be seen that the growth factor depends on the po- 
tential used. Solid line \Barreiro et al. (2000)\ dashed-line \Steinhardt et al. (1999)\ dashed-dotted line [Albrecht £r Skordis (2000)\ dotted-line 
\Brax&- Martin (1999)\ 
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Fig. 10. Predictions for Ac in the case of delayed collapse of dark energy. In all four panels, the solid line corresponds to the delayed 
collapse situation, the dashed line to the case ofT = Q both before and after turnaround. The upper left panel shows the predictions for the 
double exponential {Barreiro e t al. (2000)) , the left upper panel for the model of \Steinhardt et al. (1999)\ the lower left panel the model of 
[Albrecht £r Skordis^2000), and the lower right panel the model of \Brax £r Martin ( 1999)\ 



4. Conclusions 

In this paper we have studied the spherical collapse model in cosmologies with dark energy, provided by a scalar field 
(quintessence). Our aim was to study the model for different quintessence potentials and under different assumptions about 
the behaviour of the field during the highly non-linear regime. Our results show, that the predictions of the spherical collapse 
model depend on the form of the potential, on the initial conditions of the field and on the behaviour of the field in highly non- 
linear regions. If backreaction effects onto the field turn out to be important, these could significantly change the predictions of 
the spherical collapse model. 

We would like to point out that although the dark matter enters the highly non-linear regime first, perturbations in dark 
energy deviate only slightly from linear theory around virialisation. At the time of virialisation, the energy density of dark 
energy is always much smaller than the density of dark matter inside an overdensity. The main effect is that dark matter 
perturbations change the evolution of dark energy perturbations. The physical properties of clusters, such as the density con- 
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trast and virial radius strongly depend on clustering properties of dark energy (i.e. the parameter F), the quintessence poten- 
tial and the initial conditions for the field fluctuations. The latter is only important, however, if the field collapses together 
with the dark matter On the other hand, we have found that if the dark energy equation of state is assumed to be con- 
stant, the diff'erences between the homogeneous and inhomogeneous cases are small, at least if the equation of state w does 
not difi^er too much from w = -1. Thus, for constant equation of state the fitting formulae presented in the literature (see 
[Wang & Steinhardt (1998)" 'Weinberg & Kamionkowski (2003)^ do not change drastically even if inhomogeneities in the dark 
energy component are taken into account. 

Our results imply, that a better understanding of the behaviour of the quintessence field in highly non-Unear regions is needed 
before the predictions of the spherical collapse model can be trusted. The field can influence the dynamics of the collapse at late 
times and thus changing predictions of the turnaround, virialisation and collapse times. In order to get better predictions, the 
boundary conditions between the outer and inner metric have to be understood better. Its likely that the spherical model just gives 
qualitative but not quantitative predictions. Nevertheless, an estimate of the function T is needed, which can only be obtained 
from a general relativistic treatment. One way would be to develop a swiss cheese model in the case of a quintessential universe. 
This will be discussed in future work. 

Our calculations have also implications for models in which the scalar field couples non-minimally to dark matter. In these 
models, the scalar field perturbations and the dark matter perturbations are coupled even in the linear regime and backreaction in 
the highly non-linear regime will certainly be important. A similar conclusion might hold for models based on tachyonic fields. 

To summarise, only when the energy flux of scalar field energy density F out of a dark matter overdensity is known, the 
spherical collapse model is able to make firm predictions which can be used to make predictions for weak and strong lensing, 
the number density of clusters, etc or even for a help to use N-body simulations in dark energy models. In the case of a non- 
minimally coupled dark energy component, the predictions of the spherical collapse model will certainly be altered. 
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